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Abstract 


The connection between the quantized wave function method and the coordinate-space 
formulation is examined. The operators of the second quantization are represented in 
a sequence of configuration spaces for 1,...,2,... etc. particles. The representation 
obtained enables a simple derivation of HaRTREE’s equations with exchange. 
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Ni,N2,... 
P39,',8 
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elementary charge, 1.602 176 x 10°C 

= +, mit h = 6.626 x 10°“ Js, PLANck’s constant 
energy operator 

transformed energy operator, see Equation 26 
HAMILTON operator 

total number of particles, n particles 


operator, to which is assigned the number n of 
particles, see Equation 1 


particle 1,2,... 

summation indices 

quantized wave function 

the quantized wave function adjoint to W(x) 
wave function 

distance vector 

spin coordinate 

unitary operator 

the unitary operator adjoint to S(t) 
time 

operator for the CouLoms potential 


cartesian coordinates of all particles 


Configuration Space and Second Quantization 


The equivalence of the method of quantized wave functions with that of ordinary wave 
functions in configuration space is known in principle; however, it appears that the closer 
relationship between the two methods has not received sufficient attention. In the present 
work, the relationship between the two methods is pursued in detail. It turns out that this 
relationship is so close that one can move directly to the configuration space at any stage 
of the calculation with quantized wave functions. 


The work contains two parts. The first part has an introductory character and contains a 
derivation and compilation of known results. The transition from configuration space 
to second quantization is considered there for the case of Bose and Ferm statistics, 
whereby the uniqueness of the determination of the order of non-interchangeable factors 
is especially pointed out. The starting point for the considerations in the second part 
are the commutation relations between the quantized wave functions (wW operators). It 
is shown that these relations are satisfied by certain operators acting on a sequence of 
ordinary wave functions for 1,2,...,”,... particles. Through this the W operators are 
represented in the configuration space (more precisely: in a sequence of configuration 
spaces). Furthermore, the dependence of the w operators on time is considered and the 
form of the operator y = oe is found. Based on the obtained representation it is shown 
that the time-dependent ScHRODINGER equation for the w operators can be written as a 
sequence of ordinary SCHRODINGER equations for 1,2,...,n... particles. As a further 
application of the representation obtained, a simple derivation of HaRTREE’s equations 
with exchange is given. 
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1 Part | — Transition from configuration space to 
second quantization 


The reader who is familiar with the theory of second quantization can skip this first part and start 
reading straight away with the second part, see section 2. Relevant literature can be found in 
[1, 8,6, 5, 9] 


We denote by x, the totality of the variables of the r-th particle (e.g. the coordinates and 
the spin of the electron, x, = (%,, y;,Z,,;@,)) and consider the wave function 


W(X1, X2,..-5Xn3 f) > (1) 


which describes the entirety of n identical particles in the configuration space. It is 
convenient to move from the original variables x to new variables E with only discrete 


values 
E=E,F°)...,E™,... (2) 


through a canonical transformation; The quantities according to Equation 2 can be 
thought of as eigenvalues of an operator with a discrete spectrum. If we denote the 
corresponding eigenfunctions by 


Wx) = W(E,x) , (3) 
the transformed wave function 

(en ah ee bi ere ee 3 (4) 
is linked to the original wave function according to Equation 1 by the relation 


VC bnaMO= DS, Cie Eines li) VMimi) cl Baw) OC) 


Eij5En 
where each of the summation variables FE), F>,...,£, passes through all values of 
Equation 2. 
The Scur6piNGeR equation in configuration space is given by” 
0 
SY (ta. as) ih Se <0. 6) 


We assume that the energy operator 5 is of the following form: 


5 = Ha) + D) Ga xD - (7) 
k=1 


k<l=1 


*Here, /i denotes PLANcK’s constant divided by 27; i = 4 
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The simple sum describes the energy of the individual particles, the double sum describes 
their interaction energy. In the case of CouLoms forces 
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, € 
G(x, x’) = rey (8) 


applies. The ScHrODINGER equation for the transformed wave functions according to 
Equation 4 is obtained by introducing Equation 5 into Equation 6, and expanding the 
result according to the products 


W(E\; x1) pees »W (En; Xn) 
of the functions according to Equation 3 and setting the coefficients of the individual 


products equal to zero. There results 


> >, (EHIW) + c(E1,..., Bit, W, Bests, Bait) 
k=1 W 


+ > >) ExEIIGIWW’) 


k<l=1 WW’ 
7 c(E\, ar Rice »E-1, W, Exit, araes »Ess; W’” Eis, orace Ens t) 


O 
STs CCE pecs) Uy, 9 
i aye ( i ) (9) 
where the following designations for the matrix elements were introduced: 
(E|A|W) = | Fe9- He) wows de, (10a) 
(EE"|G|WW’) = WE; x) -W(Ex’) 
- G(x, x’): W(W; x) > W(W’; x’) dx dx’. (10b) 


The arguments 
J oe Oy aa of Pree 2 


of the wave function c in Equation 9 may be correspondingly equal to the eigenvalues 
EE"), E™, ak, EE, vate E™ . 
If we write for abbreviation 
(r|H|s) instead of (E|H|E”) , 
(rt|G|su) instead of (EO EMGIEME*) : 


C(r1,1,--+51nit) instead of c (E°), E™,...,E™; t) ; 
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then the wave equation becomes equal to 


n 
> Slbine (uate ree sine) 
k=l 


r 


+), 3 (rarilGIrs) 


rs k<l=1 
EO Ts OLS Tis ek ig Ss TIS eT) 
; 4) 
SARS SC Uis sta =O (9a) 


according to Equation 9. 


So far we have not considered the symmetry properties of the wave function, i.e. the type 
of statistics. Now the wave function (both w and c) is either symmetric (Bose Statistics) 
or antisymmetric (FERMI Statistics). In the case of the symmetric wave function, fixing 
the numbers 


N1,No,...,N;,... (11) 
is sufficient to determine the value of c(71,12,..., 13 0)°, which indicate how often the 
argument 

| ere Ae 
or 


| GAD OMe eae 
in question occurs in c. We can therefore put: 
CHE SO Ws Mh ek) (12) 


A specific sequence of numbers according to Equation 11 now corresponds to a specific 
system of values 7, 72,..., 1, regardless of the order of the latter quantities. For example, 
we have (for n = 3) 


c(4, 4,5) = c(4,5,4) = c(5,4,4) = c*(0, 0,0, 2, 1,0,0,...). 
In the normalization condition 


> e@itxcgreoP 4 (13) 


TL yveey Tn 


one can first carry out the summation over all permutations of the numbers of a fixed 
system of values 7;,72,...,/, and then over different systems of values: 


>; > le(ri,%,--..% OP =1. 


Wis.aha) Petm, 


3Jn the infinite sequence of numbers according to Equation 11 there are at most n numbers other than zero. 
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. n! 
The sum )perm, Contains > 


alr, equal terms; we therefore have 


!enyol-.. 


n! 4 
» SSS Oi =1 
ny, 7 
; 
or, if we introduce the n, as a variable according to Equation 12 


n! 2 

yy ——.. |e* (ny m,... Of = 1. (14) 
ny! * ny! Perak 

n1,N2,... 


ni 


In the normalization condition according to Equation 14 the “density function” —"~— 


can be reduced to 1 through the Ansatz* 
é* (m,m,...50) = a CN (15) 
For the new wave function f the normalization condition then follows as 
Df isms aOR =1: (16) 


Ny N25... 


In the case of Fermi Statistics, specifying the numbers n, is not initially sufficient to 
clearly determine c (7), /2,...;1), because this specification means that the quantity c is 
only determined up to the sign. But we can also retain Equation 12 and Equation 15 
for the Fermi Statistics if we introduce the additional condition that the arguments in 
c("1,12,...3;) should form a “natural” order, e.g. 


Ty <2 <13°°*+ <T. 


If the order of the arguments arises from the natural one by an even permutation, Equa- 
tion 12 holds unchanged; for an odd permutation the sign must be changed. We have 


e.g. 
c(1, 4,5) = -c(4, 1,5) = c*(1, 0,0, 1, 1,0,0,...). 


From now on we will consider the Bosr and Fermi Statistics separately. 


“For the definition of an “Ansatz’”, see e.g. https: //en.wikipedia.org/wiki/Ansatz 
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1.1 Bose Statistics 


In the case of Bose Statistics, several identical arguments can occur in the wave function 


C(71,12,...51), e.g. 
C = C(U,U,U,V,V,W,...)- 


In the first sum of the expression according to Equation 9a, functions can therefore occur 
that only differ in the order of the arguments, namely there are n, terms in which the 
argument r takes the place of u, n, terms with r instead of v, etc. If we combine the same 
terms, we get the expression 


> (u|A|r) n, - c(r,u, U,V, Vv, W,...) 


+ >» (v\A|r) ny - c(u,u,u,r,Vv,w,...) 
| i 


feee 


for the first sum in Equation 9a. If we introduce n,; as a variable here according to 
Equation 12, we get 


> UHI) ty C8 (ote Lye tet 1) 


+) WHI) ny c* (ty Leet td.) 


+.--- 


or more simply 
DD (PII tp c* (....mp— Lem td), (17) 
Dp r 


where the index p now can pass through all values (and not just the values p = u,v, w,...), 
since the superfluous terms disappear because of the factor n,. For r = p, 


CP (castip = lien ly) 


simply means 


Analogously one can also transform the second sum of the expression in Equation 9a. If 
we take into account the number of identical terms, we get: 


1 
>; {wuts . aru -(n, — 1)-c(r, 8,U, Vv, Vv, W,...) 


r,S 


+ (uv|G\rs) -nyny, > c(r, uu, 8, V,W,...) 


1 
+ (vvIGIrs) » Sy» (my — 1) -clunar si. , 
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and if we introduce the quantities c* (n;,,...): 


1 
>, {autre = 1) C8 oeaty=2yeeeat Dyecett lee) 


rs 


+(uv|G|rs) -nny-c* (...,n, —1,...,m—1,...,0, +1,...,n5+1,...) 


1 
+ (wIGIrs)- Sn (My = De" Go. ty = 2,225 teams + 1...) 


ve} 


or more simply 


; . > > (pq|GIrs) - ny - (n, - Spa) 


D.q 18 


2O (esata Teac Ieccodn et Ty sen eh Lyac) 5 (18) 


Here again the summation indices p and gq can pass through all values without excep- 
tion (and not just p,g = u,v,w,...). The factor !/2 must be present in all terms, since, 
for example, both the combination p = u and gq = v as well as p = v and q = u oc- 
curs in Equation 18. The meaning of the term in Equation 18 for the case that two or 
more of the numbers p, q, r, s coincide probably does not require any special explanation. 


Using Equation 17 and Equation 18, the wave equation according to Equation 9a can be 
written as: 
DD, (IHN) + tp c* (...5 1p - 1,...,n,+ 1.52) 
Dp r 


+ yy > (pqlGIrs) - np - (1, 7 pa) 
Pq 


rs (19) 


ger (ae tipTedsuslin= Leeson Tous eos) 


. 0 
~ ih ac" (Mi Ma... = 0. 


To proceed further, it is convenient to introduce the operator U,, which transforms a 
function 


T Gist Mist) 
into 
f(y, N2,--. Ny + Lge) 


according to 


O64 F (ghee he) Sy (istic oe (20) 
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The matrix of U, and its adjoint U * are of the form 


(21) 


oo Oo Oo 
oo OF 
Co oF Oo 
oF Oo O&O 
S 
ll 
Co oF Oo 
oF Oo OO 
jet ED! CE RO 
oo Oo oO 


with respect to the variable n,. Consequently, the adjoint operator U; transforms the 


function f (11, 72,...,;,...) into f (m1,N2,...,n,- 1...) ifn, # 0 and into 0 for n, = 0. 
So 
Ulf (1,M2, Ms) =F my... — L(t #0) ox 
=0 (n, = 0) 
From the definition of U, it follows 
U'-U,=1. (23a) 
On the other hand, U, - Ui # 1, namely 
0 0 0 O 
0 1 0 O 
O25 =O). A (23b) 
0 0 0 1 


So the operator U,, is not unitary. 
Furthermore, for p # r U, and U; are interchangeable with U, and Uy. 


Using the operators U,, the functions c* appearing in Equation 19 can be written in the 
following form: 


ral Ore acetone rae Mee 0 A oug @iey (rear norma! 
Oa hltg = lyapen ge Vancdtth Leena lines) 


—yiyi * 
=U,U UU, -¢ (itisttp, crip ian niu aMpis Mics 


The order of the factors (U* to the left of U) follows clearly from the definition of c* for 
p =r inconjunction with Equation 23a and Equation 23b. These expressions apply to 
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arbitrary (even coincident) values of p, g, r, s. If one introduces them into Equation 19, 
one gets 


oS > (pI Ir) -n,USU, - c* (ni, 1, ...) 


rs 


P q 
1 
ie » )) (palGirs) - np « (mq = bpq) « UZULU.U, + c* (msm, «) 


: (6) 
~ ih Fe" m,...30 = 0. (19a) 


Here we still have to express c* (m1, 2,...) according to Equation 15 by f (,7,...). 
The operator n for the total number of particles, and consequently also n!, apparently 
commutes with the products Uj, -U, and Uj, -U -U,.- Us; further we have 


1 
VWny!-nol-... 


1 1 
VN. Nge:... 


n!} 
ny!-ng!.... 


tion 19a is therefore equal to 


| n! [n!-ny!- 

; try 1 Qhteee 
ml-ny!-... npU,U, n! f (m1,M,...) 
= yn, ULU,- yn,- f (usm,...) . 


Analogously, using Equation 24a and Equation 24b and the relation 
(14 - 5p) Ut = Utn, 


-U,+ ni! -ng!-...= Vn, +1-U, =U,- Vn, (24a) 


The one with 


multiplied term npU; U,c* (nj, N2,...) of the first sum in Equa- 


for a term of the second sum in Equation 19a, we obtain the expression 


n! nyl-no!-. 
SS ae = -ut 5 eet eRe, 
Vata 5nq) ULU UU, > 4] = Ff Gjstinyccs) 
= yn,-U)- yng: Ul -U,- Vns-U,+ Vn; f (n2,-..) « 
If one introduces these expressions into Equation 19a, one obtains for f (7, 72,...;31) 


the wave equation 
Of _ 


Hf (m1,Mo,+..30) -1h- 2 =0, (25) 


where H denotes the transformed energy operator 


H =) (plHIr)- yapULU, vi, 
pr 


1 
+ 5° Dy (palGlrs) - yitpUy WigU, Us VigU, Vir (26) 


pars 
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The operators U, and n, only appear here in the combination 


b, = U,- ro 
(27) 
Bi sain, U] 


If one introduces Equation 27 into Equation 26, one gets the following expression for H: 


1 
H = 575}, (pir) by + 5») 6}; (palGlrs) bsby (28) 
pr 


paqrs 


As can be seen from the definition according to Equation 20 and Equation 22 of U, and 
U;, the just introduced operators b, satisfy the relations 


bi -b, =n, 
b,- bi =n, +1, ”) 


and furthermore, for r # s the b, and bi are interchangeable with b, and bi, then one has 
the well-known interchange relations 


b, : bi = bi : b, = Ors, (30a) 
b,-b,-—b,+b, =9. (30b) 


If one now uses b, to form the quantized wave function 


Wx) =D by Wel) (la) 


with its adjoint 


W(x) = DPE VCO, (31b) 
the energy operator H can be represented in the following form: 
H = ih Wi (x) - A(x) - Wx) dx 
+ ; { Wi (x) W(x’) GO, x’) WO’) W(x) dx da’. (32) 


The commutation relations for the quantized wave functions (w operators) follow easily 
from Equation 30a and Equation 30b taking into account the equation 


DV) WAX) = d(x). 


One obtains the expression 


Wx’) W(x) — W(x) WX’) = 5(x - x’), (33a) 
W(x’) W(x) — Wx) - WO’) = 0. (33b) 
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1.2 Fermi Statistics 


We return to the wave equation Equation 9a. We think of the numbers 7,72, 73,... 
ordered in the natural order 


re <%<173<+9+ <M, (34) 
so that we have according to Equation 12 
CU Pisesss he) Re ON ee). (35) 
In the natural order, the number 7; is at position 
kK=nt+innt+-:++n,.- (36) 


In the k-th term of the first sum in Equation 9a r; is replaced by r, so that the arguments 
in c are there in the order 


VYSV 2,0 eo Tk-15 Ty Vhtly ++ eo Tn (*) 
which is no longer a natural one: Here r is at position k, while it should be at position 
Fe EP , , 
ko=n,+n,+--->+n,. 


(The deleted quantities are the new values of the n,, which correspond to the arguments 
according to Equation * in c.) Therefore it is 


kak 
CVs is TRB = CY HO Gast ae the) 


It is useful at this point to introduce the operators @, and a} by setting 


caignece isis) tong = 0; 

Ue Gcccti sae, ONE ks (37a) 
0 for n, = 1, 
fe ray, 

POP Giese ley (37b) 
Ff(m,...,n,-1,...) for n,=1. 


From this definition it follows that for r # s the operators a, and a} are interchangeable 
with a, and a (since they act on different variables) , while for r = s the equations 


Ee ys 
(38) 


Q;- a =1-n, 
hold. Furthermore, one can easily prove the equation 


a,:(1 - 2n,) = —(1 - 2n,)- a@,. (39) 
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Using the operators a, one can write 
* * 
OU aes Ies eles) = al a,-c (N,N...) . 


The order of the factors a, and a, is clearly determined here, because for r = 7, and 
n,, = 1 the factor of c* on the right must reduce to 1. We have 


(- 1) = (-1)0 tH 


and since for n = 0 and n = | the quantity (—1)” agrees with (1 — 2n), we can also write 


for it: 
Tk 
(-1)' = | | (1 -2n,) = vn, 
p=l 
where 
vs =| | (1 -2n)) (40) 
p=l 


is the WicneR’ sign function. Analogously, 


(-D* =, 


r? 


where v;, is formed with the numbers n/.. So we have 


, 


3 = 7 x 
Cy 150355 Tits Tete SD) = Vn VO Gre © (1, Md,.05.)s 


and because of 


vaia,=a'a,v 
eel) ne es | ee ae 


we can write for it: 
lg \ t * 
OW hi eesT iste Teas 252) = Veg OVO Op Nas <3) 


The first sum in Equation 9a is therefore equal to 
n 
> pS (rlHIr) + V0} ary, + C* (M1, M2,-..) « 
r k=l 


When summing over k, the index r; runs through the values 7% = 11, /2,...,1n. Instead, 
one can let 7; pass through all values, because the superfluous terms disappear due to the 
property of the operator a’. So we get the expression 


» os (pIHIr) vpata,y, + c* (ny, M2...) (41) 
Dp a 
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for the sum under consideration. We now want to transform the second sum in Equa- 
tion 9a. Above all, it is about determining the sign in the equation 


BOA foccg les Tielg ees Tote Ss igs has) 
* 
=O Gig tin = laccstty = weit tgs A) 


a ok (er, een ee 


We first transfer the argument r, which is in the k-th position, to the first position in c; 
This gives c the factor —(—1)* = —v,, and we have 


Chips tay Mite hissy MAGS Piles tag Te) 


SS U5 Figs FT eies ss ST Sai TO : 
If 7; > rz, s remained at position 
b=ny tng t+ sng 


for 7 < rg, however, s has moved one place to the right and is now at position / + 1. If 
we now transfer s to the second position, we get 


—VyVrC (V, 8,115,125 +++) for r; > rz, 
CO... Meets Vo Vht ly ee 5 Meds Ss Vla1s +++) = 
+Vn Vine (T, S,11,12,...) for TL <x. 


On the other hand, if we now denote the natural positions of r and s by k’ and /’ according 
to 


=n +n ,+---+n,, 
Van, tnyt-e--+n,, 
we get 
natural order 
EL eT eto a Sipe) CO Noyes Bek aes oe ae LOL es 
t+ViVi.C(1,5,11,1%,.-.) fors<r 


by a completely analogous consideration. Together with the previous equation, this gives 


COs Fets hs oT Ss Tis ee) 


+V 7, Vr ViVC* (n\, Rosie ) in case I, 


s~N OSS 


—Vi,VnVVC* (n\, ee ) in case II, 


where cases I and II are characterized by the inequalities 


Tr) > Vr ands >r 
or case I, 


r<r,ands <r 
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rT, > rz ands <r 
or case II. 
rT) <x ands >r 


Replacing the arguments 1,,7,... with n},nj,... in c* is caused by the operator 


Past 
Q},Q},AsQ;: 
* , , * 
6" (js) = ai a) aa, -c* (N,N, ...) 


One can be convinced that the order of the factors a‘ and a (as far as it is important) is 
chosen correctly by considering the special cases r = r,, s = 7, and r = r;, s = rz. If we 
also take equation 

VV a} a} as, = a! a) Aa,ViVs 


into account, we have 


OG ii Thats BH Peigoea Tis Tig) 


$V Vr, OAV sV/C¥ (n1,N,...) incase I, 


Vy, Vrj Qh, Oy, sQ,V sV/-C* (nj,n2,...) incase II. 
Now it follows from Equation 39 and from the definition according to Equation 40 of v, 


@-Vs = vVs-a@,forr>s, 
(42) 


Q-' Vs = —-v,:a@,forr<s. 
Therefore, in case I we have either 


Ay’ Vs = Vs + A, 


and forr>sandrn>r, 


at the same time or 


y+ Vs = —Vs* Oy 
and forr< sandy <r 
t ed t 
On, Vir = Vr, * Ur, 


at the same time. So, in case I, the operator acting on c* (11,2, ...) is equal to 
Ae t a 
VO, Vin; OL sVsQrVy « 
But this operator also has the same sign in case I, because we then have either 


Wy * Vs = Vs * Oy 


and forr > sandr <r 


t = f 
On, Vy = Vr ry 
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or 
@, + Vs = —Vy° 
and forr<sandrn>r. 
a}, "Vr = Vr," a}, 
So we always have 


a foe at * 
Cos Hes Feels ev TA Sia ees) = Vy, VO, OV Oe CNT, Moy es4) 


We now have to introduce this expression into the second sum in Equation 9a. This sum 
will be equal to 


n 
»S », (reF GIs) + Vp Q}, Vil AsV sQVy + CX (M1, N2,--.) » 


rs k<l=1 


If we drop the restriction k < / here, the sum doubles and we have to add the factor !/2. 
We then get 


rs 


1 
5 s >; (pq|G\rs) - Vp) VGA hQsV QV CPN Mig Nara) (43) 
Pq 


Here p and q initially only pass through the values 71, 72,...,7, (where p # q). But one 
can let them run through all values without exception if one takes into account that the 
superfluous terms disappear. 


Substituting Equation 35, Equation 41 and Equation 43 into Equation 9a gives the wave 
equation for the wave function c* (1, 12,...; ). However, in the case of Fermi Statistics, 
this wave function differs from the wave function f (7;,72,...) according to Equation 15 
by only one factor (namely vn!), which is interchangeable with the individual terms of 
the energy operator. Therefore, the wave equation for f (7, 72,...) has the same form as 
for c* (n1,N2,...), namely 


Hf (mm,..s)-in SD =0, 


where the energy operator H after Equation 41 and Equation 43 has the following form: 


1 
H = >, (pIA|r) - Vpah ary, + 5 » (pq|GIrs) - Vp VGA hQsVsQrVp (44) 


pr pqrs 


In the energy operator H, the operators a, and v, only appear in the combination 


ay = Ayr' Vr, 
i t (45) 
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namely 
1 
a t Tat 
H = » a, : (plH |r) - a, + a y aa, (pglG|rs) - asa, . (44a) 


paqrs 


As can be easily proven with the help of Equation 38 and Equation 42, the “quantized 
amplitudes” a, satisfy the equations 


Bap = Tp, (46) 
a,-a) =(1—-n,) 


and the commutation relations 


apa t+a'-a,=O6rs; (47a) 
a,*d,+da,:da,=O0. (47b) 

If one uses the “amplitudes” a, to form the quantized wave functions 
W(x) =D) ay ol), (48a) 
Wix)= di ai U0), (48b) 

which satisfy the commutation relations 

Wx’) WI) + W(x) WO’) = (x= 2), (49a) 
W(X!) Wx) + WO) WO’) = 0, (49b) 


then one can write the energy operator in the form 
H = { W(x) A(x) - W(x) dx 
1 
a a Wi (x) WO’) Gx, x) WO’) - WO) dx dx’, (50) 


just as in the case of Bose Statistics. 


The transition from the quantized amplitudes a, (or b, in the case of Bose Statistics) to 
the quantized wave functions w(x) represents a unitary canonical transformation of the 
variables of the one particle (Transition from the E” according to Equation 2 to the x). 
The a, (or the b,) can be viewed as quantized wave functions just as well as w(x), and 
the formule represented using the a, (or b,) (like e.g. the commutation relations or the 
expression for the energy operator) are equivalent in content to those written with the 


W(x). 


It should also be noted that all other operators that can be represented in the configuration 
space can also be transformed according to the pattern of the energy operator and 
represented using the quantized wave functions. Just as with the energy operator, the 
order of the non-commutative factors results in an unequivocally clear manner. 
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2 Part Il — Representation of the / operators in 
configuration space 


In the formule of the second quantization, the total number n of particles no longer 
appears explicitly; the formule apply to any or even an indefinite n. The number n can 
be assigned the operator 


n= [ ¥'Q)-worde, (1) 


which has the eigenvalues n = 0,1,2,... 


All operators can be divided into two classes with regard to their behavior towards the 
operator m: the first class includes those that are commutative with m, and the second 
that includes the non-commutative operators. An independent definition of the operators 
of the first class is given in the concurrent work by P. Jorpan, [7]. It is shown there 
that all such operators can be constructed by the coordinate-space-method in a form 
which is completely identical for both statistics; The exact same formule even apply to 
the mathematically possible but physically inconceivable solutions to the many-body 
problem. It then follows that commutative rules apply to all of these operators, which are 
independent of the statistics. 


Here we want to deal with the representation of the general operators in the configuration 
space that are not commutative with m; Above all, it is the representation of the operator 
w(x). Of course, the results can also be applied to the operators that are commutative 
with m, since these can be expressed as w(x) and w*(x). 


In order to capture both types of statistics uniformly, we write the commutation relations 
for the quantized wave function in the form 


Wx’) W(x) — W(x) WO’) = (x - 2), (2a) 
W(x’) Wx) — €- Wx) - Wr’) = 0, (2b) 
where for the Boss Statistics ¢ = +1 and for the Fermi Statistics e« = —1 is to be set. 


From the definition according to Equation | of the operator m and from the commutation 
relations according to Equation 2 it follows that 


m-y—w-(m—1)=0 (3) 


for both types of statistics. We choose a representation for w(x) in which m has a 
diagonal form. If we denote the matrix elements of w(x) in this representation by (n|W|n’), 
it follows from Equation 3 the “selection rule” 


(n—n' +1)-(nlp|n’) = 0, (3a) 
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which states that only matrix elements of the form (|y|n + 1) are different from zero. 
The matrix W(x) is therefore of the form 


0 OWL 0 0 
soe 0 OO ClW2) 90 (4) 


0 0 0 (2|y|3) 


A single matrix element (n — 1|y|n) in Equation 4 can be viewed as an operator acting on 


a function of n variables® x1, 22,...,X, and transforming this function into one of (n — 1) 
variables x), .x2,...,X,-1; and the parameter x. So, the operator w(x) acts on a sequence 
of functions 
const. 
(x1) 
W(x1, X2) (5) 
W(xX1, X2, X3) 


of 0,1,2,3,... variables and transforms them into an analog sequence; namely, the 
functions according to Equation 5 can be understood as ordinary SCHRODINGER wave 
functions in configuration space®. We want to say that W(x), x2,...,X,) 1s the wave 
function in the “n-th subspace”. We further want to show that the commutation relations 
according to Equation 2 are satisfied by the Ansatz 


(1 — 1 y(x)|r0) We, x2, - a) = Wn W(X, 1, X04. 6s Xnt) (6) 


if one defines the adjoint operator w'(x) starting from w(x) correctly. According to 
Equation 4 the matrix for (x) is of the form 


0 0 0 
Wi = (110) 0 0: . (4a) 
0 (21) 0 


where (nly*(x)In - 1) is the operator adjoint to (n — 1|W(x)|n), which transforms a func- 
tion of (n — 1) variables x,,x2,...,X,-, into one of n variables x;,x2,...,x, and the 
parameter x. It should be noted that the operator (nly*(x)In — 1) must not change the 
symmetry property of the wave function: it must convert a symmetric function into a 
symmetric one and convert an antisymmetric into an antisymmetric. We want to find the 


Each “Variable” x, is actually the set of variables, such as x,, y-, Z-, 0, Which describe the r-th particle. 
Such function sequences were first considered by L. Lanpau and R. Peterts, [10]. 
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operator (nly*(x)In — 1) by determining its kernel; we can do this if we form the kernel 
of (n — 1|y(x)|n) and then move on to the adjoint kernel. 


Since the wave function is either symmetric or antisymmetric, one can also write 


(n — 1|y(x)|n) W (x1, XQ, 2025 ay) 


IWAN Noses og Ane PE WO Ns oe A) 


1 
I a { 
$e bE) (XL D5 Xn ¥)] (6a) 
instead of Equation 6. The kernel of the operator defined by Equation 6a is 


(n = 1; x1, Xo, see Xn W(x) |n3 1, &2, coe onl 
1 
= —-{6(€; —x)-d0(&—X1)-...° 6 (E, — Xp 

A toe )-6(&-m) (En Xn-1) 

teeth! 6(E — x1) + 20.6 (Gea — Xe) O(E - x) 

6 Spa — X%) +. O(En — Xn-1) 

Hertel 5 (E, — x1)... 6 Ent — nt) 6 En — Oh. (7) 
The kernel of the adjoint operator (n|y(x)|n — 1) is obtained by replacing €,, &,...,&, by 
X1,%X2,...,X, and x1, X2,...,Xp-1 by €1,&2,...,&,—-1 in Equation 7. The result of applying 
the operator (nly*(x)In — 1) to the wave function wW (x1, X2,..., Xn-1) is therefore equal to 


(nly (x)|n = 1) : W (x1, XQ5-065 Xn) 


1 
— — - {5 (x — X)° W(X, .X3,...,X,) E+ 0 (X. — X) W(X], X3,..-5 Xp) 


vn 
ee ONES A) UA Ape ho esis hpyio wane tp) 
eee 8 y= 2) Ws, %at)} - (8) 


The operator (nly*(x)In 7 1) defined by this equation satisfies the requirement that it 
leaves the symmetry property of the wave function unchanged; the transition from Equa- 
tion 6 to Equation 6a was carried out precisely for this purpose. 


After w(x) and w*(x) are defined, we can proceed to the proof of the commutation 
relations according to Equation 2a and Equation 2b. We have to form the operators 
Wi (x) - W(x’) and W(x’) - wi(x). These operators are interchangeable with m, so have 
diagonal form with respect to n. We have 


Ag O O 
wig-wo@yat® Ar 0 (9a) 
0 O A 
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and 
Bo O O 
wxy-wayafO Br Ot (6b) 
0 O B, 
where A, and B,, denote the operators 
Ay = (nl'(x) - WO )In) = (ab Coin — 1) «(2 = LO)In) (10a) 
and 
By = (nlWx’)- WOIn) = (IWO)In + 1) - (2 + LY In) (10b) 


which act in the n-th subspace. By first applying Equation 6 and then Equation 8 we find 


An W(X], X2,..., Xn) = 0 (x1 — x) W(X, %2,...5X%n) te 
groan a Cae emeee Mire seem 
$25 (ty — 2) WO Met), ” 


or, if we take into account the symmetry property of the wave function, 
(nlyh" (x) WOX’)In) «Wr x2, Xn) 
= 6(4) — x) W(W',%0,..., Xn) te +O (MR - 2X) 


Ue 6 Mg 
te +6 (Xp, — X) W(X, %2,---5 Xp, X') « (11a) 


If one first applies Equation 8 and then Equation 6 and replaces n by n + 1, one finds: 
By W (x1, X2,-+-5Xn) = (nlp(x’) yw (x)In) WH Xi, Hog 3k ey Xp) 
= 6(x' — x) W(x, X0,..-, Xn) FEO (%1 — X) WH, X0,. Xn) Hoo 


+ BF <6 Ole — 0) WR, Mig se os Kes Khas so In) HO 
+ €"-6(X, — x) W(X, X1,---,Xp-1) - (12) 


Comparing Equation 11 and Equation 12 shows that 


Ba + W201, X25 06+, Xn) — & An W(X, X2,.--%n) =O(xX—- 2’) 
CW Os Nog cg ahs (13) 


According to Equation 9a and Equation 9b this means that the commutation relation 
Wx") W(x) — € W(X) WO’) = 6(x- x’) (2a) 


exists, where the identity matrix (with respect to n) must be added to the right. 
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Proving Equation 2b is even easier. According to Equation 6, the operator w(x) transforms 
the sequence of functions according to Equation 5 into 


W(x) 
V2 - W(x, x1) 
V3 + W(x, 1,3) 
Le. 
const. W(x) 
W(x) V2 + W(x, x1) 
W(x) =) Wr.) P= 4 VB WO )> - (14) 
W(X1,X2,%3)] fete 


If one applies the operator w(x’) to Equation 14, one gets 


const. v2-1- W(x, x’) 
W(x) V3--2- W(x, x’, x1) 
W(X’) Waxy =4 WX.) 7 = V4-3- W(X, x’, X1,X2)¢ - (15a) 
W(x1,X2y%3)| | te ee ee ees 


By exchanging x with x’ one gets 


const. V2 -1- W(x’, x) 
W(x1) V3 -2- W(x’, x, x1) 
W(x) WO’) =) Worx) 6 =) V4-3- WO, xx1, 32) ¢ - (15b) 
WX1,X2,X3)] | tee eee 


Now the expressions according to Equation 15a and Equation 15b are either equal 
(€ = +1, symmetric functions) or oppositely equal (e = —1, antisymmetric functions), 
which proves Equation 2b. 


With the help of the formule obtained, all operators of the second quantization can be 
constructed using the coordinate-space-method. Those of them that are not commutative 
with m act on sequences of functions of the kind according to Equation 5 and cannot 
be represented in a space of a certain number of dimensions. For the operators that can 
be commutated with m, it is sufficient to consider the diagonal element of the matrix 
with respect to n, which can be viewed as an operator in the n-th subspace, i.e. in the 
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configuration space for a fixed number nv of particles. For example, the energy operator 
according to Equation 50 is commutative with respect to m, and its construction in the 
configuration space leads back to the usual ScHRODINGER energy operator for n particles. 
Here we want to look at a few more examples of the operators that can be commuted 
with respect to m. 


The operator w(x) - W(x) of the particle density has the representation 


(kb (x) - W(x)In) Wr, x2, «5 Xn) 
= {5 (xy — x) $5 (2) $+ + 6 — )} WO rn) - (16) 


in the n-th subspace. This formula is a special case of Equation 11, which is obtained by 
setting x’ = x in Equation 11 and using the relation 


6 (x% — x) > f(x) = 6 (x, — x) f(x) 


valid for every continuous function f(x). 


If we multiply Equation 16 by the volume element dx and integrate over a subvolume V, 
we can conclude that the operator 


my = { Wi (x) W(x) dx (17) 
V 
in the n-th subspace has the following representation: 
(n| my |r) - Wr (X1, X2,.- 2, Xn) = My (X1, X2,-6 6 Xn) W(X, X02, Xn) (18) 
The value of the function ny, (x1, X2,...,X,), which is multiplied in Equation 18, is equal 
to the number of arguments x), x2,..., x, that belong in the subvolume V. The operator 
(n| my |n) therefore has integer eigenvalues nj, = 0, 1,2,...,n, as required. 


As a further example, we consider the operator for the CouLoms potential 
(x) - W(x D4 
vey = f SE ¥ (19) 
lr — v’| 
To form its matrix element (| V(x)|n), we replace x in Equation 16 with x’, multiply by 


e 
and integrate over x’. We obtain 


jr—v’] | 


n 2 
CHIVCOIN) “(18a = DS Weir) 20) 
k=1 


ee 
The operator V(r) therefore means a “multiplication by )17_, — 


iy) in the n-th subspace. 


26 


2 Part II — Representation of the w operators in configuration space 


2.1 Time dependence of the & operators and the quantized wave 
equation 


As is well known, the change in the state of a physical system over time is expressed 
either in the time dependence of the wave functions or in that of the operators. According 
to Dirac, [2], we refer to the representation of the operators that corresponds to the 
time-dependent wave functions as the ScHRODINGER and the one corresponding to the 
time-dependent matrices as the HEISENBERG depiction. 


Let y be the wave function of the system and S(f) be the unitary operator’, which converts 
the initial value W(-,0) of the wave function into its value W(-, f) at time t. (The variables 
of the system are indicated here with a dot.) We have 


W(-,t) = S(t) - W(-,0). (21) 
If one differentiates this equation according to time and replace Y(-, 0) by 
WC0) = S10) We, (21a) 
one gets 
OW. : ; 
9p 50° WO.0 =SOSO We). (22) 


We denote the Hermitian operator i S'S‘ as 
: : 1 
iS()S'(t) =-iS(HS"(t) = - KH, (23) 


because SS* = 1. # is then the Hammron operator of the system. If we denote by L 
an operator in the SCHRODINGER representation, and by L’(t) the same operator in the 
HEISENBERG representation, we have 


L'(t) = S"(OLS (0). (24) 
From Equation 24 and Equation 21 there results 
L'OWC,0) = S"(NLW, 1). (25) 
The time derivative of this expression is e to 
* Oy,0) =< (SLM) (26) 
a oO | 


On the left is the operator —~ in the HEISENBERG representation. If we denote the same 


operator in the SCHRODINGER repieschtaion with wee 


and Equation 25 


, we have, analogous to Equation 24 


dL’(t) 
dt 


7see https: //en.wikipedia. org/wiki/Unitary_operator 


dL 
= Si(t)- Zo (27) 
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and 
dL’(t) 


dt 
Comparing Equation 26 and Equation 28 yields 


dL 
W-,0) = S*(t)- Pres (28) 


dE i 28: LISOLNK 
God = 80-5 (SOLU) . (29) 


or, if one performs the differentiations, 


N= SOS OMI Dw a 
Gp 8 OS OLD EHC. D) (30) 


Note that according to Equation 21 (or according to Equation 21a) the operator S'(f) (or 
S(t) causes the temporal continuation of the wave function around time ¢ in the positive 
(or in the negative) sense, the content of Equation 29 can be formulated as follows. 


In the ScHRODINGER representation one gets the result of applying the operator ve to the 
wave functions W-, f) if one carries out the following operations: 


1. Application of the operator L. 

2. Time continuation in negative time direction around time ¢. 
3. Differentiate according to time. 

4. Time continuation in positive time direction around time ¢. 


This formulation has the advantage that the HAMILTon operator is not used explicitly. 


This rule can now be applied to the determination of the operator ow which appears in 
the quantized wave equation. The operator L in this case is the quantized wave function 
W(x, t) and W(-, t) is the function sequence according to Equation 5. We limit ourselves to 
the case where the number of particles does not change over time and thereby exclude the 
consideration of photons. To obtain more descriptive formule, we consider the quantized 
SCHRODINGER equation 


a 
[H°Cx) + Vad] WO) = in- a (31) 


where H°(x) denotes the ordinary ScurOpINGER operator of the one-body problem: 


h2 2 2 2 
hanes (2 Pe a 


+ =] + U(x, y, Z) (32) 


2m \ae2” dy” a2 


and V(x) = V(2) is the operator for the CouLoms potential defined by Equation 19. 
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In our case, the operator S(t) in the configuration space simply causes the temporal con- 
tinuation of the individual wave functions W(x, X2,...,%Xn3t) of the sequence according 
to Equation 5 


const. const. 
W(x1; 0) W(x1; t) 
S(t)4 W(x1,%230) p={4 Wx, x30 7. (33) 


W(X, X2, X33 0) W(X), X2, X33 1) 


So the operator S(t) has the diagonal form 


Son Or 1G) 6 
0 Sr 0 0 

ROS 0 On. 10 Eps (34) 
0 0 0 S; 


where S,, = S,(¢) is an operator that continues the wave function in time in the n-th 
subspace: 
S nt) > Ws (x1, X2,....%05 0) = Wa, X2,.-- nD (35) 


We also have according to Equation 23 
—if-S,OS}OQ = H (x1, %2,...,%n) (36) 


where # (x1, X2,...,X,) denotes the HamILTon operator in the n-th subspace. The opera- 
tor —ih-SS' is also a diagonal operator with Equation 36 as diagonal elements. 


op 


We now form the operator w(x; t) = F 


. We have according to Equation 29 or Equation 30 


const. 
W (x13 0) 
Wt) 4 W(x, x25 0) 


W (X1, X2, X33 0) 


w (x; 0) W(x; 1) 
V2- W(x, x13) " V2 - W(x, x51) 
V3 W(x, x1, 4050) V3 + W(x, x1, X25 2) 


= §(S*(1) (37) 
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and if we consider Equation 36 


const. 
. . W (x13 0) 
h- : 
ie W (X1, X25 0) 
0 W(x; 0) 
— V2 - H(x1)W (x, x15) _ V2 - W(x, x13) a8) 


| a SG Gael VNB Gea 


For the operator on the left side of the quantized wave equation according to Equation 31, 
we get, if we take Equation 20 into account, 


const. 
Hees W(x13 1) 
[H°(x) w+ Vox) | Fee 
W(x; t) 
V2 - W(x, x13) 
= H?® V % i 
| ae o) V3 + W(x, x1, X25 0) 
H°(x)W(x; t) 
V3: |H%() + | Wx n138 
= 2 y (39) 
V3 - [ae us nea — W(X, X1, X25 £) 


By equating Equation 38 and Equation 39 we get (after dividing by V2, V3, etc.) a 
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chain of equations 


Hx) Worst) = in EOS (40,) 
0 e . fe) 
A(x) + + FE(x1)| W(x, x50 = 1h: —W(x, x15), (402) 
lr = r|| Ot 
e e 
[He a Hs.) 
pee je =| 
: 6) 
W(x, x1,%3 0 =ih- ae X1,X230), (403) 
n e 
09 + > o + P(01, Xa... =) 
|e - tl 
. (6) 
WG Nips cass das) = 2 As Hp Ms Ms Hs «+ Ani t) (40n41) 


From Equation 40 one can see that 
HEX) = A(x); 


is satisfied; Equation 40, then returns 


2 
€ 


IE (x, x1) = H°(x) + H° (x1) + etc. 


lr — ry| 


In general, the (n + 1)-th equation provides a recursion relation between the SCHRODINGER 
operator for n and for (n + 1) particles, namely 


n 2 
e 
FO (Xy Xs X25 +005 %n) =H (x) + Oe Eis Bip 2 REY (41) 
oS |r - tl 
If one now expresses H# (x1, X2,..., Xn) directly in terms of H’ ° one gets for the HAMILTON 


operator of the n-body problem the common ScHRODINGER expression 


n 


FO AX HRS Ra) = SH (xy) + oy 
k=l 


k>l=1 


2 
€ 


. (42) 
Ite — Tl 


The quantized wave equation, Equation 31, therefore breaks down in the configuration 
space into a sequence of ordinary SCHRODINGER equations according to 
., OW 
FO (Hy Hasse Hn) WA aes Met) HS An- a (43) 
This example shows that computing with quantized wave functions allows an immediate 
transition to the usual configuration space at every stage. 
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2.2 Derivation of Hartree’s equations using the method of second 
quantization 


As a simple application of the results obtained, we want to derive the HARTREE equations, 
[4], supplemented by taking the exchange into account. 


As is well known, the equations for the eigenfunctions of the energy operator (as well as 
the HARTREE equations) can be derived from the variational principle 


OW =0, (44) 


where W denotes the energy of the atom in the considered stationary state. It is therefore 
sufficient to find the expression for the energy. But now W is equal to the diagonal 
element 


W = (Wn| H|Wn) (45) 


of the matrix for the quantized energy operator H (see Equation 50 of the first part): 


ae FOeaeie 
2 af W(x) WO’) WO’) - WO) WX) 4 


(46) 
peep 


To determine W, we need to calculate the matrix elements of the operators under the 
integral sign. 


We have 


(Wahp'(x) - H(x) - WO)|Wn) = A(x’) - (Wahb') - Wor)|Wn) (47) 


for x’ =x. 


To calculate the matrix element of the operator in the first integral, it is sufficient to 
determine the quantity 


o(x,x) = (Walp (x) YO)|Wn) . (48) 


We have already found the expression for the operator y"(x) - w(x’) in the n-th subspace 
(see Equation 11a); With the help of the eigenfunction Ww (x1, X2,...,Xn) of the energy 
operator in the n-th subspace, which belongs to the eigenvalue W, we get from this if we 
take into account the symmetry property of the wave function, 


o(x, x’) wine foe Pye eats te) aso) da diy (49) 
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To derive the Hartree equations, we must replace the wave functions in this exact 
expression with an approximate expression in determinant form according to 


1 
Ww (%1, X2,---5Xn) = val - |ei(xx)]| (i,k = 15253515050) (50) 
where the y,(x) are assumed to be orthogonal and normalized: 
ie Pix)» r(x) dx = Gig . (51) 
We then get 
O(x,x') = DB) gil’), (49a) 
i=1 


and consequently, according to Equation 47 and Equation 48, 
(Walp (x) - Hx) - WoW) = D1 Gx) - Hx) - gil). (52) 
i=l 


We now calculate the matrix element of the operator in the double integral in Equation 46. 
Using Equation 14, Equation 16 and Equation 8 one easily finds that this operator has 
the following meaning in the n-th subspace: 


(alex) WO!) WO’) WOdIn) «WO, X25. Xn) 
= DO te = 2) 6 (1 = 2) «rs a5 en) (53) 


k i=l 
(k#l) 


Consequently, its matrix element is equal to 


(Walp (x)- W'O")- WO’) -WOo|Wn) 


= ne(n= tye fof ow (aa! aa...) 


After inserting the determinant expression according to Equation 50 for ww, the approxi- 
mate formula is obtained from Equation 54 (cf. [3]) 


(Waly) WO’) Wor) WOoIWn) = ofa.) - oC’, - loxP. 5) 


If one now introduces Equation 52 and Equation 55 into Equation 46, one gets for the 
matrix element of H, i.e. for the energy W the expression 


dx3:--dx,. (54) 


W = ip DF) HO) #1) dx 


2 e. [fee x) O(x’, x) — los, x’)P diac’. (56) 


lr — 1’| 
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This expression only differs from the one given in our cited work, see [4] (formula (93)), 
in that here we consider the spin coordinate included in the variable x and are therefore 
entitled to operate with purely antisymmetric wave functions. 


Leningrad, Physical-Technical Institute. 
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